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Abstract 

By using k-factorization, we derive resummation formulas for the non-abelian qq 
contributions to both heavy flavour production by gluon fusion, and to the next-to- 
leading BFKL kernel. By combining this result with previous ones by Fadin et al. 
on the virtual terms, we also compute in closed form the complete qq contribution 
to the gluon anomalous dimension in the Qo-scheme. We find that qq resummation 
effects are important for heavy flavour production, but are instead small in the 
anomalous dimension eigenvalues, because of a cancellation between abelian and 
non abelian contributions. 
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1 Introduction 



Hard processes in the small-x regime s > Q 2 > A 2 play an increasingly important role 
at present colliders ||] and are characterized by the fact that the effective QCD coupling 
constant a s {Q 2 ) log(l/x) is sizeable and the corresponding perturbative contributions 
need to be resummed. 

Such resummations are performed, at leading log level, on the basis of the BFKL 
equation || for the anomalous dimensions, and of the k-factorization formulae [3-6] for 
the coefficient functions. Particular attention has been devoted in the literature to heavy 
quark production processes [Q, [(J and to the (abelian) light quark contributions to the qg 
and qq entries of the anomalous dimension matrix§. Furthermore, an impressive program 
of evaluation of next-to-leading (NL) kernels of the BFKL equation is under way [7-10]. 

In this note, we focus our attention on the (non-abelian) qq contribution to the k- 
factorization program, which is relevant in two ways: firstly, as a hard subprocess in 
Drell-Yan heavy-flavour production and secondly, in the massless quark limit, as an im- 
portant NL contribution to the gluon anomalous dimension. We shall then derive QCD 
resummation formulas in both cases and in particular we shall obtain in closed form all NL 
terms of the form a s Nf(a s CA^og(l/x)) n , contributing to the gluon anomalous dimension 



in the so-called Q - sc heme [11, [12|. 



The evaluation of such contributions is particularly important for structure functions. 
Since the gluon couples to the proton only through qq states, the NL anomalous dimension 
7 99 (determined essentially by the abelian qq kernel) contributes to scaling violations at 
the same level as 7 9r It is then important to perform a complete calculation of NL qq 
contributions (in particular, the non-abelian one), in order to check the pattern described 
above, and to reduce the factorization scheme dependence of the theoretical estimates of 



scaling violations at HERA |L3, [L4 . 
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2 QQ Hadroproduction 



Let us start considering the cross-section for the hadroproduction of a (heavy) quark- 
antiquark pair, which at high energies is dominated by the gluon fusion process in Fig. 1. 
This contribution may be expressed in the k-factorized form ||: 

MV= f ^^d 2 k 1 d 2 k 2 a(k 1 M,M 2 ,z 1 z 2 s)^ 1 \z 1 ,k 1 )J r{2) (z 2 ,k 2 ), (2.1) 
J Z\ z 2 

where T denotes the unintegrated gluon density in the hadron and a the high-energy 
projection of the (Regge) gluon fusion process g(ki)g(k 2 ) — > QQ. 

Since the evolution of the structure functions is simpler in the space of z and k moments 
it is useful to express the high energy factorization formula (2.1) in terms of the double- 
Mellin transformed structure functions 

and hard cross-section coefficient function [] 

so that it takes the form: 

ljr9 ( M 2 \ fa s \ fh+^.dj-td^duj ( s \ w ( M 2 \ 71+72 ,,,, > , , , n 

Ua {-H-)i,J J ^{jp) (of) ^ 1, ^)^(7 2 )^(0,,7 2 ).(2.4) 

For large enough M 2 , the 7 integrals will be dominated by the BFKL anomalous 
dimension, 71 ~ 72 ~ 7x(a s /a;). This means that the high-energy resummation effects on 
the cross section we are considering are embodied in the 71, 72-dependence of the hard 
coefficient function ^(71,72), which is precisely what we wish to compute. 

-"^This function is related by the change of normalization h — -^7172-ff to the /i-functions of Refs 
ill 
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The squared matrix element for the process under consideration was computed in Ref. 
(see also Ref. ||), and contains two terms, with colour factors Cp and Ca, that we 
shall call respectively the abelian and the non abelian contributions. 

The abelian contribution H a is known |J and provides the quark anomalous dimension 
||. Let us then concentrate on the non abelian one H na . 

First, we rewrite the squared matrix element of Ref. || in terms of the Sudakov 
parametrization for the exchanged gluons' momenta 

k? ~ z lP 1 + k£, e 2 ~ z 2 p% + k£, (2.5) 



and for the momentum transfer 



A» = z 1 x 1 p f t - z 2 x 2 p^ + A", (2.6) 



where p%, p 2 denote the (light-like) momenta of the incoming hadrons. 
By using explicit expressions for the invariants 



8 = {h + k 2 ) 2 , i = A 2 , u = (h -k 2 - A) 2 , v = s + (ki + k 2 



(2.7) 



the non abelian squared matrix element of Ref. 0] can be rewritten in the form 



Ana _2„ 

A = Ti a 



* h - — 7. (1 - Xi - x 2 ) H h 

{M 2 -t){M 2 -u) s\M 2 -u M 2 -t) y 1 us 



+ kfk| U + 2 (1 - - x 2 )- 

k 2 (l - x 2 ) + k 2 (l - a*) - ki ■ k 2 + A ■ (k 2 - k 2 ) 



1 x x x 2 v v k 2 (l-x 2 )+k^(l-xi)-k 1 -k 2 + A-(k 2 -k 2 ) 
2"M^-2 (1 - Xl - X2)+ i 



(2- 
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By then integrating over the phase space and performing the moments of Eq. (2.3) 
we obtain the non abelian coefficient function 5™ a (7i, 72) as follows 

^(71, 72) = ^ /-T^kf — ^^(^(1 - xOi/ - (k! - A) - M ) x 

/ k 2 \ 71 / k 2 \ 72 

x5(x 1 (l-x 2 )^-(k 2 + A) 2 -M 2 ) f-^J f-^J A™(ki,k 2 ,A,Zi,^,M 2 ) (2.9) 

This expression can be evaluated analytically (at least in the limit u = 0, relevant at high 
energies) with a careful choice of the order of the integrations. We find it convenient to 
start eliminating the variable v by integration of one mass-shell delta function, and then 
computing the A-integral with the aid of a Feynman parametrization of denominators. 
The remaining integrals, first over the transverse momenta ki, k2 and then over the 
Feynman parameters and longitudinal momentum fractions, can be evaluated in terms of 
Gamma and Beta functions, and give, after some algebra, the result 



#" a (7i,72) 

_ C A Pis 
2 7T 



r(i - 71 - 72)1X71)1X72) 



5(1-7!, 2 -705(1 -72, 2 -72) 



(3-2 7l )(3 - 2 72 ) 



^ 5(1 - 7l ,l- 7l )5(l- 72, 1-72) 
(l + (l_ Tl )(l- 72 ))j + 



+5(71, 1 - 71)5(72, 1 - 72) 
+2 



r(2- 7 i-7 2 ) r(3- 7 i-7 2 ) 



r(4 - 271 - 2 72 ) r(6 - 271 - 2 72 

5(71, 1 - 71 - 72)^(72, 1 - 7i - 72)5(3 - 71 - 72, 3 - 71 - 72) 



+ 



(1 - 71 - 72) 



(2.10) 



The details of the calculation are reported elsewhere . 



Let us note the triple pole singularity at 71 + 72 = 1 of the last term of Eq. (2.10), 
which is related to the collinear-singular behaviour of the partonic cross section in the 
massless limit. In fact, if M 2 <C (ki + k2) 2 = q 2 <C k 2 ~ k 2 , the hard cross section 
approaches the singular limit 



^ =0 (ki,k 2 ,M 2 



_ a s f 1 xi(l-xi) / xi(l-xi)q r 

~-a s — / dxt - log 1 +■ 

71 Jo q 2 \ 



a. 



M 2 
a s C A 



~a s — 


log 


q 2 


5" 


/ D7T 




M 2 


~ 3 



71 



(2.11) 



which, after taking the k^-moments, yields precisely a triple pole. The factor in front of 
the logarithm is related to the g — > qq GLAP splitting function |3]. 

Such singular behaviour was already discussed in Ref. H and yields an interesting 
enhancement of the cross section when 7x and 72 approach the BFKL saturating value 
7 = 1/2. Here we further notice that, after combination with virtual terms in a BFKL 
kernel, the behaviour (2.11) is responsible for the qq contribution to the running coupling 
constant, as we shall see in the following. 

Indeed, the same cross section b occurs, in the massless limit, as a light quark contri- 
bution to the BFKL kernel, given by the inverse of Eq. (2.3), i.e., by 

' ; ^M,M>) 1 t~ ^H^JKY 11 m\ (2,2) 



a s M 2 "" u W " J M 2 Ji-ioo (2tt) 2 uv ' x ' iaj \M 2 J \M 2 _ 

After subtraction of the singular contribution of Eq. (2.11), the eigenvalue of the finite 
part of the kernel is directly obtained by the massless limit of Eq. (2.12), which singles 
out the residue at the (single) pole at 71 + 72 = 1, as follows 



fCi(7)= lim (tf na (7,72)-#r(7,72))(l-7-72) 

72— ►1—7 



7T J 2 



r 2 (7)r 2 (i-7)a-(2 + 37(i-7)/ (2 - 7)r(1 + 7) 



(2.13) 



^3 v ,v ,y r(4-2 7 )r(2 + 27) / 
for each quark flavour. 

Note that the various terms of Eq. (2.13) show spurious double poles at 7 = and 
7 = 1, which would signal the presence of a collinear singularity in the t-channel. They 
cancel however in the sum, as they should because the real emission non-abelian diagrams 
do not contribute to the t-channel singularity. 



3 Next-to-leading gluon anomalous dimension 

Let us now come to the important issue of NL contributions to the flavour singlet anoma- 
lous dimension matrix 7 a h(a s ) (a, b = q,g). As is known, the gluon entries start at leading 
level 

Igg = 1L (^jj +a s lNL (^jj + •-, 7 OT = ^7l(1 + 0(a s )), (3.1) 
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while the quark entries |J start at NL level and, in the Qo-scheme are given by 



% g = N fl lH 2 ( lL ), l gq = ^(l q9 -^). (3.2) 

Here if 2 is the k-factorization vertex measured by the _F 2 structure function and, similarly 
to Eq. (2.13), is given by the abelian qq contribution as || 

H 2 ( 7 ) = Cj: 1 lim fff( 7j72 )(l-7-7 2 ) = 

72^1 — C 

« s r 2 ( 7 )r 2 (i -7)r(2-7)r(i + 7 ) , , „ . . 

T r(4-2o0r(2 + 2 7 ) 7 ' (2 + 37(1 - 7) ). (3.3) 
The NL terms of 7 99 occur in the larger of the anomalous dimension eigenvalues 

. C F C F N f a s 
7+ ^ Igg + 7^7 9S , 7- ^ (3-4) 

and require a full discussion of the BFKL kernel 

K u (a„ k 1? k 2 ) = (K (k u k 2 ) + ^(ki, k 2 )) . (3.5) 

UJ 

Here i^o is the leading kernel, whose 7 -dependent eigenvalue 

Xo(7) = 2^(l)-V(7)-^(l"7) (3-6) 
determines the anomalous dimension 7 i(o7 s /u;) 



™ \ 2 



1 = -Xo(7l), 7l = - + 2C(3)(^) +... . (3.7) 



Furthermore, both Kq and the NL kernel K\ are scale-invariant, because the non- 
invariant NL part has been factored out in Eq. (3.5) inside the running coupling constant 



at the scale ~k\ [Tj| [17 1 



a s (kl) = a s (fi 2 ) - ba 2 s (n 2 ) log ^ + /i = O(Q). (3.8) 
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We shall compute the NL contributions (which are factorization scheme dependent JTT|, 



18|), in the so-called Qo _sc heme of Ref. ||12|| , further discussed in a forthcoming paper 



fl5| . In this scheme we define the quark sea density as in the physical DIS scheme |T9 
i.e. by the structure function F2, similarly to the parton model. 

Instead, the gluon density is defined by the procedure of k-factorization itself, which 
provides a gauge invariant off shell continuation, so that the initial gluon virtuality can 
be fixed at Qq. Therefore, gluon states carry a high-energy polarization which in general 
does not match the collinear definition at higher orders. For this reason, the expression 
of 7 gg in Eq. (3.2) does not carry the R-factors of the MS type schemes [|ll], ||, and a 



leading similarity transformation [12| is needed to relate the former scheme to the latter. 



The key advantage of the Qo-scheme is that the anomalous dimension 7+ of Eq. (3.4) 
is computed from the eigenvalue of the scale invariant kernel Kq + a s K\, directly in four 
dimensions, as follows: 



ex 

1 = — (Xo(7+) + asXi(7+)) , d s = d s (fi 2 ). (3.9) 

UJ 

By expanding and using Eq. (3.7), we obtain the NL contribution to 7+ 

,NLf®s\ __ „ Xi(jl(u s /u)) _n L . C F 



Vo;/ xo{lL{a s /u)) C A 

which contains, in particular, the qq contribution. 

In order to compute Xi a -> ^ ne non-abelian part of xVi we need to combine the singular 
real emission contribution in Eq. (2.11) with the virtual corrections, and then to add the 
resulting finite term to the nonsingular real emission part in Eq. (2.13). 

The virtual qq contribution to the one gluon emission kernel was computed in Ref. 
and, after azimuthal averaging, takes the form 

a s N f ( 1 , M 2 1 , k?\ , . , 

-irU l0 ^-kFk| l0g k| )■ (i=k 1+ k 2) (3.11) 
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in which fi 2 = 0(Q 2 ) is the factorization scale, and M is the (small) quark mass. By 
combining Eq. (3.11) with the singular qq emission part in Eq. (2.11) the M 2 -dependence 
cancels out and we get 



a s N f 

67T 



log ^ + log v2 



q- 



5 
3 



k? 



k? 



df log k 2 



(3.12) 



A similar calculation can be performed |15j directly at M = and D = 4 + 2e 
dimensions, to yield the ^-dependent kernel 



a s N f 

671 



i r 2 (l + e) /q 2 
eT(l + 2e) \/i 2 



8 q 2 k 2 - k 2 10g k 2 



(3.13) 



In order to regularize the q 2 = singularity in Eq. (3.12) and (3.13), one should 
combine the above results with the purely virtual term, computed in Ref. ||10||. This can 



be done in dimensional regularization by means of Eq. (3.13), and the result [T5| can 
be read off directly in four dimensions as follows: the logarithmic term at scale ki in 
Eq. (3.12) is absorbed in the running coupling factor of Eq. (3.5) and the remaining 
(regularized) collinear part becomes 



a s (K[ v) + K? 



a s N f 
6tt 





-!) 








q 2 



k? 



R 



kf 



k| l0g k 2 



(3.14) 



where we have introduced the notation 



/(k 1; q) \ R = /(k 1; q) - 5 2 (q) rf 2 q/(k 1 , q) 



(3.15) 



and q 2 > A 2 is understood in all q-integrations. It is now straightforward to diagonalize 
the kernel (3.14) and to combine it with the finite real emission part in Eq. (2.13) to get 
the complete non-abelian eigenvalue expressed in terms of Xo(j) of Eq. (3.6) and of VKt) 



N f a s 

67T 



71 



sin 7T7 



1 — 3 cos(7T7) 



1 + §7(1-7) 



(3-2 7 )(l-2 7 )(l + 2 7 ) / 



(3.16) 
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together with its abelian counterpart of Eq. (3.3): 

(a) C F a a N F ( 71 V l + f 7 (l-7) 

°' Xl = Un^J C ° s(7r7) (3-2 7 )(l-2 7 )(l + 2 7 ) - (3 - 17) 

Note that both eigenvalues are symmetrical for 7 — > 1 — 7 , apart from the first term in 
square brackets of Eq. (3.16) (~ x'o) which comes from the fact that we have factorized 
the running coupling effects in the upper mass k^. 

In order to extract from Eqs. (3.16) and (3.17) the anomalous dimension eigenvalue 
7+ we have to use Eq. (3.10), i.e., we have to divide it by x'o Since x' ( 7 )~ — 1/ 7 2 + 0( 7 ) 
for small 7 , the one-loop and two-loop results in the DIS scheme are read off from the 
small 7 behaviour of Eqs. (3.16), (3.17), that is 

(a), N ~ C F a s N f ( 1 13 , \ ( n0 ) ~ a s Nf( 1 23 A 

asXi (t) ki^ — + — + 0(1) , a s xi + - ( 3 - 18 ) 

Ca 37T y 7 z 6 7 / 67r y 7 Z 6 7 ) 

In fact, by using Eq. (3.10) we check directly the known |21J low-order expressions in 
the DIS scheme 



Furthermore, Eqs. (3.16), (3.17) provide an all-order resummation in the Qo-scheme. 
Since 7(?g is known from Eq. (3.2), and — Xo7 2 = 1 + 0( 7 3 ), we find that x? is n °t 
perfectly cancelled by (Cp/CAjjqg, and therefore 7 ^ L takes also contributions of type 
Nf Jq^* s (ff) from three loops on. Finally, the terms Nf 7t ° !s are resummed by x'" a , 

which is typically of negative sign, as it appears already at low order in Eq. (3.19). 

The total qq eigenvalue xVil) — Xi(7)+Xi a (7) is plotted against 7 in Fig. 2. It 
is apparent that a cancellation between abelian and non abelian terms is at work, and 
that resummation effects are small, in their sum, even for 7 values which approach the 
saturation limit 7 = 1/2. Only for 7 — > 1 is the abelian term dominant, because of its 
stronger 7 = 1 singularity. 

This remark is made more significant by the observation that the scheme dependence 
of 7+ is a simple one. For instance, the transformation from the DIS-Qo scheme to the 
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DIS-MS scheme, involves a scale change by a factor R('yi J (a s (t))) Ullj under which 
7 + changes by running coupling effects only, i.e., 

7+ = 7+ - |- (3.20) 

We tentatively conclude that qq resummation effects in scaling violations are small 
in a combination proportional to T q + (Ca/Cf)^ which isolates the 7+ eigenvalue, and 
is accessible experimentally through proper combinations of F2 and -Fl[|T|, fL8|| , when the 
latter will be measured. It is important to check this conclusion against the purely gluonic 
NL terms, which are still under investigation M, [|. 

A final comment is needed for the running coupling effects (to which the qq channel 
contributes) which remain an important NL feature of the BFKL equation. We have 
decided to factorize the a s dependence at the upper scale kf Jl7j in Eq. (3.5). It can 
be shown [[L5[] that, in the regime k 2 ^> Qq 3> A 2 , with d s (Ql)/u; < 1, this dependence 



is equivalent to a renormalization-group evolution due to the anomalous dimension 7+ of 
Eq. (3.9), apart from a calculable coefficient factor, common to quark the gluon in the 
Qo-scheme, given by the expression 



7l (jjj V~Xo(7l) 



(3.21) 



Since the latter differs from 1 only from 3 loops on, the choice we have made is good 
at perturbative level. The factor (3.21) can be reinterpreted in terms of an additional 
contribution to the gluon anomalous dimension, along the lines of Eq. (3.20). 

To sum up, in this paper we have resummed in closed form the non abelian qq con- 
tributions to both heavy flavour production and to the small-x anomalous dimensions. 
While non abelian resummation effects are quite important in the former case, they nearly 
cancel with the abelian ones in the latter, due to virtual contributions in the anomalous 
dimension eigenvalue. 

The most significant NL effects (apart from the running coupling) seem to remain the 
rather large scaling violations due to 7 99 flT3"| , which somehow represent a leakage from 
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the gluon density to the quark sea, which unfortunately depends on how the gluon is 
defined, even in the class of DIS-schemes. It is therefore important to measure a quantity 
independent of F2, for instance Ft,, in order to disentangle the gluon density from the 
quark sea and/or to look for combinations with smaller scaling violations, as is suggested 
by the result on 7+ of the present work. 
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